1. IntroductionThe iron-based superconductors have been investigated in depth both experimentally and theoretically.[1–7] Most of the iron pnictide superconductors have strong thermal fluctuations compared to most of the other superconductors (except cuprate superconductors) due to their large Ginzburg–Landau (GL) parameter κ, large critical temperature Tc, large anisotropy parameter γ, etc.[8] For different structures of iron pnictides, such as (1111), (111), (11), and (122) families, different dimensionality behaviors in the superconducting fluctuation region have been claimed.[9–16] However, the discussions for the dimensionality are still inconclusive. It was claimed that the superconducting fluctuation of LiFeAs followed the 2D behavior above
in Ref. [11], but it was questioned in Ref. [17].
To describe the iron-based layered superconductors, we use the Lawrence–Doniach (LD) model, in which the original bulk expression is divided into two-dimensional (2D) superconducting layers with a Josephson coupling between them. The LD model describes well the experimental data of cuprate superconductors, for example, the data of the fluctuation conductivity and diamagnetisation.[18,19] Due to different temperatures and different magnetic fields, there exist 3D region, 2D region, and crossover region, which are shown in the vortex liquid state of the H–T phase diagram in Fig. 1. Region I (below Tc) is the 3D region. In region I, the coherence length along the c-axis is larger than the layer distance. Taking into account the 3D limit, where the layer distance can be taken to zero limit in the theoretical formula, the theoretical conductivity curves are almost the same as the original ones. Region III is the 2D region, where the coherence length is quite small compared to the layer distance. Region II is the crossover region between region I and region III, where the coherence length is on the same order with the layer distance. The larger anisotropy parameter, which leads to strong thermal fluctuation, may enhance the 2D and the crossover regions. Hence the strong anisotropy of the iron-based superconductors can enlarge both regions II and III in the H–T phase diagram compared to the less anisotropic superconductors.
Both phenomenologic and microscopic multi-band theories have been developed to describe the experimental observation of the multi-band electronic structure in most of the iron pnictides families, such as (1111), (11), (111), and (122) families.[20–22] However, the single-band model can be used approximately near Tc.[23] We use the single-band LD model in this paper to study the fluctuation conductivity of vortex liquid near Tc. The often used method describing the fluctuation properties of superconductivity is Gaussian fluctuation theory (GFT),[13,15,24,25] in which the quartic term of the GL free energy is omitted. The GFT method can only describe region III in Fig. 1. The self-consistent approximation[18,19] was developed and was used to explain the various experimental data, like conductivity and diamagnetisation of vortex liquid. In the self-consistent Gaussian theory (SCGT), the quartic term is taken into account in the Hartree–Fock approximation. The SCGT can describe all three regions in Fig. 1.
This paper is organized as follows. In Section 2, the calculation of the superconducting fluctuation conductivity based on the time-dependent Ginzburg–Landau (TDGL) function with the linear response approach under the self-consistent approximation is briefly introduced. In Section 3, we present the theoretical fitting of the experimental conductivity data and the discussion. We summarize and conclude in Section 4.
2. Fluctuation conductivity based on the TDGL modelThe mean-field GL free energy of the LD model is
![](cpb_25_9_097401/cpb160960eqn1.gif)
where
d′ is the distance between layers labeled by
n,
m* is the effective Cooper pair mass in the
ab plane, and
mc is the effective Cooper pair mass along the
c axis. The anisotropy parameter
γ is denoted as the ratio of the effective masses along the two different orientations,
![](cpb_25_9_097401/cpb160960ieqn3.gif)
.
![](cpb_25_9_097401/cpb160960ieqn4.gif)
, where
![](cpb_25_9_097401/cpb160960ieqn5.gif)
, in which
![](cpb_25_9_097401/cpb160960ieqn6.gif)
is the mean-field critical temperature. The covariant derivative is defined by
D = ∇ + i2
π/
ϕ0A. The superfluxon is
ϕ0 =
hc/2
e with electric charge
e = |
e|. The ratio of the coherence length
ξ2 =
ħ2/(2m*
αTc) and the penetration depth
λ2 =
c2m*
b′/(16
πe2αTc) describes the Ginzburg Landau parameter
κ =
λ/
ξ, which is very large for iron pnictides. We will consider a constant external magnetic field. Due to large
κ, the diamagnetization is
B −
H ∝
H/
κ2 →
B ≃
H. Hence we can choose
A = (−
By,0), where
B =
H, and assume that the magnetic field is orientated along the crystallographic
c axis. The supercurrent in the original unit is given by the GL theory
![](cpb_25_9_097401/cpb160960eqn2.gif)
2.1. TDGL modelTo describe the transport properties of the layered iron-based superconductors, the TDGL equation will be used, i.e.,
![](cpb_25_9_097401/cpb160960eqn3.gif)
where
Dτ = (
∂τ − i(2
e/
ħ)
Φ) is the covariant time derivative, in which the scalar electric potential is
Φ = −
E′
y. The electric field is therefore assumed to orientate along the
y axis. The inverse diffusion constant
γ′/2 is assumed to be real, as we can neglect the tiny imaginary part of
γ′, which has to be included when the Hall conductivity and Seebeck effect are considered.
[19] ζn is the thermal noise,
[26] which satisfies
![](cpb_25_9_097401/cpb160960eqn4.gif)
where ⟨⋯⟩
th is the thermal average. We substitute Eq. (
1) into Eq. (
3) and obtain
![](cpb_25_9_097401/cpb160960eqn5.gif)
where |
ψn(
r,
τ)|
2ψn(
r,
τ) will be substituted by 2⟨|
ψn(
r,
τ)|
2⟩
ψn(
r,
τ) using the Hartree–Fock approximation, and ⟨⋯⟩ means that we take thermal and space-time average if not specified.
The physical units will be used to make the model and equations dimensionless. The coherence length ξ will be used as the unit of length and the upper critical field Hc2 ≡ ħc/2eξ2 as the magnetic field unit. The dimensionless order parameter is
. The TDGL equation will take the dimensionless form
![](cpb_25_9_097401/cpb160960eqn6.gif)
where
ε = −
ah + 2⟨|
ψn(
r,
τ)|
2⟩ and
ah = (1 −
tmf −
b)/2. The covariant derivative is
D = (
∂/
∂x − i
by)
i +
∂/
∂yj, in which the dimensionless magnetic field is
b =
B/
Hc2. The dimensionless covariant time derivative is
Dτ = (
∂τ + i
Ey), in which the electric field
E =
E′/
EGL is scaled with the unit
EGL =
Hc2ξ/
cτGL.
τGL is a typical “relaxation” time in the superconducting phase,
τGL =
γ′
ξ2/2. The layer distance
d (
d =
d′/
ξc) is rescaled with the coherence length in the
c axis,
![](cpb_25_9_097401/cpb160960ieqn8.gif)
. The dimensionless Gaussian white noise takes the form
![](cpb_25_9_097401/cpb160960eqn7.gif)
where
ωmf is the fluctuation strength,
![](cpb_25_9_097401/cpb160960ieqn9.gif)
, and
![](cpb_25_9_097401/cpb160960ieqn10.gif)
is the customary Ginzburg number,
![](cpb_25_9_097401/cpb160960ieqn11.gif)
.
The ⟨|ψn(r,τ)|2⟩ can be calculated by solving Eq. (6), and ε can be obtained self-consistently by the gap equation[18,19]
![](cpb_25_9_097401/cpb160960eqn8.gif)
The dimensionless supercurrent density takes the form
![](cpb_25_9_097401/cpb160960eqn9.gif)
with the unit of the current density,
jGL =
Hc2c/2
πκ2ξ. The conductivity will be given in units of
![](cpb_25_9_097401/cpb160960ieqn12.gif)
.
2.2. Calculation of the conductivityThe linear response theory is used to calculate the superconducting fluctuation conductivity σs. The total electrical conductivity is σ = σn + σs, where σn is the normal conductivity.
To solve the linearized TDGL Eq. (6), we use the Fourier transform
![](cpb_25_9_097401/cpb160960eqn10.gif)
The TDGL equation becomes
![](cpb_25_9_097401/cpb160960eqn11.gif)
where
![](cpb_25_9_097401/cpb160960eqn12.gif)
Then the TDGL equation is rewritten as
![](cpb_25_9_097401/cpb160960eqn13.gif)
where
ψkz(
r,
τ) can be expended to the linear order of the electric field, i.e.,
![](cpb_25_9_097401/cpb160960eqn14.gif)
and
![](cpb_25_9_097401/cpb160960eqn15.gif)
in which
j(
kz) is denoted as
![](cpb_25_9_097401/cpb160960ieqn13.gif)
. The zero order and the first order equations are solved as follows:
![](cpb_25_9_097401/cpb160960eqn16.gif)
![](cpb_25_9_097401/cpb160960eqn17.gif)
The supercurrent in the
y direction is obtained in the Fourier form
![](cpb_25_9_097401/cpb160960eqn18.gif)
We substitute Eqs. (
16) and (
17) into the supercurrent Eq. (
18), and only keep the terms with zero order and first order of
E. Then the supercurrent takes the form
![](cpb_25_9_097401/cpb160960eqn19.gif)
The first term in the above equation is the supercurrent of the equilibrium state, which equals to zero, and we only need calculate the second term. In order to calculate
![](cpb_25_9_097401/cpb160960ieqn14.gif)
, the following eigenfunction of
![](cpb_25_9_097401/cpb160960ieqn15.gif)
will be used:
[24]
![](cpb_25_9_097401/cpb160960eqn20.gif)
where
Lx is the sample length along the
x direction,
t* is the time period,
k,
Ω, and
l are discrete quantum numbers, and
![](cpb_25_9_097401/cpb160960eqn21.gif)
Equation (
18) can be solved by the insertion of complete set ∑
l|
φl⟩⟨
φl|, in which all the Landau levels are added. Here we take
l = (
l,
k,
Ω) for simplicity. The supercurrent density becomes
[24]
![](cpb_25_9_097401/cpb160960eqn22.gif)
where the factor
![](cpb_25_9_097401/cpb160960ieqn16.gif)
is used to take space–time average of the supercurrent, and the Dirac bra–ket
![](cpb_25_9_097401/cpb160960ieqn17.gif)
of the quantum mechanics matrix element is used for facilitating the calculation. The matrix elements in Eq. (
22) need to be solved first. Since
φl is the eigenfunction of operator
![](cpb_25_9_097401/cpb160960ieqn18.gif)
, the matrix elements
![](cpb_25_9_097401/cpb160960ieqn19.gif)
,
![](cpb_25_9_097401/cpb160960ieqn20.gif)
, and
![](cpb_25_9_097401/cpb160960ieqn21.gif)
can be calculated easily with Eq. (
20) and quantum mechanics of Landau levels.
Discrete summations can be transferred to integration
. The supercurrent therefore takes the form
![](cpb_25_9_097401/cpb160960eqn23.gif)
and the integration of
Ω is
![](cpb_25_9_097401/cpb160960eqn24.gif)
The supercurrent becomes
![](cpb_25_9_097401/cpb160960eqn25.gif)
in which cutoff (
Nf) of the Landau levels is imposed. The final result is
jy =
σsE, and
![](cpb_25_9_097401/cpb160960eqn26.gif)
Here the ultra-violet cutoff is defined as
Λ = (
Nf + 1)
b and the cutoff
Λ is the depairing energy of Cooper pairs. Compared to the inter- and intra-Landau-level excitation energy scales in dimensionless units, i.e.,
b and
ε,
Λ is sufficiently large near
Tc, so the expression of the supercurrent can be weakly cutoff independent.
σs is consistent with the solution of the Green function method, when cutoff
Λ approaches infinity.
[19] 2.3. RenormalizationThe self-consistent equation ε is solved by substituting Eqs. (14), (16), and (17) into Eq. (8)
![](cpb_25_9_097401/cpb160960eqn27.gif)
The ⟨⋯⟩ can be calculated following the similar steps in the last subsection, and the result is
![](cpb_25_9_097401/cpb160960eqn28.gif)
which was derived in Ref. [
18] where the thermodynamic properties were studied based on the GL model.
The real critical temperature Tc can be defined at ε = 0 when the magnetic field is zero. By taking the limit b→ 0, the above equation gives
![](cpb_25_9_097401/cpb160960eqn29.gif)
where
![](cpb_25_9_097401/cpb160960ieqn23.gif)
,
![](cpb_25_9_097401/cpb160960ieqn24.gif)
, and
![](cpb_25_9_097401/cpb160960ueqn1.gif)
Substituting Eq. (
29) into Eq. (
28), we obtain the gap equation
![](cpb_25_9_097401/cpb160960eqn30.gif)
where
![](cpb_25_9_097401/cpb160960eqn31.gif)
and
![](cpb_25_9_097401/cpb160960ieqn25.gif)
with
t =
T/
Tc. It is easy to prove
ωt =
ωmftmf. The
ε in Eq. (
30) can be solved self-consistently, and it shall be positive definite due to the positive requirement of the variational free energy.
[18] 2.4. 3D limit and renormalizationBy taking the 3D limit of Eq. (26),
![](cpb_25_9_097401/cpb160960eqn32.gif)
The self-consistent equation in the 3D limit is therefore
![](cpb_25_9_097401/cpb160960eqn33.gif)
where
U0(
d → 0) and
Uε(
d → 0) are
![](cpb_25_9_097401/cpb160960eqn34.gif)
by taking the limit of
d → 0.
3. Experimental data fittingThe superconducting fluctuation conductivities of both the F-doped LaOFeAs compound (Tc = 19.8 K, 1111 family)[27] and the single crystal BaFe2−xNixAs2 (Tc = 20 K, 122 family)[25] are fitted in Fig. 2. The fitting parameters are listed in Table 1.
Table 1.
Table 1.
![](https://cpb.iphy.ac.cn/html_resources/images/table-icon.gif) Table 1. The fitting parameters for LaOFeAs and BaFe2−xNixAs2. .
Material |
Gi |
γ |
Hc2 |
Tc |
ξc |
d′ |
LaO(F)FeAs |
0.0015 |
7.64 |
50 T |
20 K |
3.35 Å |
8.717 Å |
BaFe2−xNixAs2 |
0.00002 |
2.076 |
61.65 T |
19.8 K |
11.1 Å |
6.36 Å |
| Table 1. The fitting parameters for LaOFeAs and BaFe2−xNixAs2. . |
3.1. LaOFeAsThe layer distance d′ of LaOFeAs used in the calculation is d′ = 8.717 Å. The expression of relaxation time γ′ in BCS,
, is derived for BCS superconductors. For high-Tc superconductors, we denote
, where constant
may not equal to one and will be determined by data fitting. The conductivity σ = σs + σn, reciprocal to the original experimental resistivity data, is fitted using Eq. (26). σn = 2.86 mΩ−1 ·cm−1 according to the experimental data.[27] The cutoff Λ is set to be 0.3, which is approximately in the same order as used in Ref. [25]. The absolute value of Λ is not sensitive to the calculation in regions I and II (see Fig.1) due to ε, b ≪ Λ.
The fitting parameters are chosen as
,
, the anisotropy parameter γ = 7.64, and the Ginzburg–Landau parameter κ = 73.36, which are basically consistent with Refs. [28] and [29]. Based on the fitting parameters, Ginzburg number
and the fluctuation strength parameter
is in the same order of magnitude with that of cuprate YBaCuO,[19] which suggests that they may have the same order of fluctuation strength. The dimensionless layer distance d = 2.6 for LaOFeAs demonstrates a relatively layered behavior.[9] Furthermore, the mean-field critical temperature
K, calculated with Eq. (29), is close to the onset critical temperature measured in experiment
, which is relatively large in iron pnictide superconductors. Therefore, the crossover region and the 2D region, i.e., regions II and III in Fig. 1, are larger than those of the other iron pnictide families, like BaFe2−xNixAs2, which will be discussed later. From Fig. 2(a), we can see that this material has a large temperature transition region, which is corresponding to the 2D–3D crossover region in Fig. 1. So LaOFeAs shows a layered characteristic in a relatively wider range of temperatures and magnetic fields compared to BaFe2−xNixAs2.
3.2. BaFe2−xNixAs2For BaFe2−xNixAs2, the layer distance is d′ = 6.36 Å. The normal conductivity σn = 0.027 × 107 Ω−1 ·m−1 is taken to be constant by ignoring the tiny dependence of σn on the magnetic field as shown in Ref. [25].
The fitting parameters are
, κ = 52.9, and Hc2 = 61.65 T, which are consistent with Refs. [25] and [30]. Thus, the Ginzburg number Gi = 0.00002 is rather small compared to that of LaOFeAs, indicating a relatively weak fluctuation effect for BaFe2−xNixAs2. The dimensionless layer distance d = 0.57 indicates that BaFe2−xNixAs2 has a 3D characteristic. Furthermore, when we take 3D limit d → 0 in the numerical calculation, it is hard to find any difference from the previous result, suggesting a 3D behavior for BaFe2−xNixAs2. Besides, the mean-field critical temperature
calculated from Eq. (29) is
, which is much smaller than that of LaOFeAs. Therefore, the crossover region and the 2D region are much smaller than those of LaOFeAs, and it shows that BaFe2−xNixAs2 is dominated by the 3D behavior (region I in Fig. 1).
3.3. DiscussionIn order to clarify the dimensionality of the ion-based superconductors, we need to compare the coherence length and the layer distance in the c axis ξc(T,B). ξc(T,B) is dependent on the magnetic fields and temperatures, which is not easy to obtain from the theoretical equations. Dimensionless coherence length in the c axis ξc(T,B)/ξc is approximately in the same order of magnitude as
, in which ε is calculated numerically from Eq. (30) with the parameters listed in Table 1. Comparing ε and 1/d2, we find that the iron pnictides can be classified by their dimensionality, as shown in Fig. 1 for three regions, i.e., region I is of the 3D feature, region III is of the 2D feature, and region II is the crossover region.
The information about the coupling between the nearest layers can be described by the Josephson coupling parameter,
. Γ is of the order of
, where α is of the order of 1.[31] When εd2 ≫ 1, the coupling between the layers is weak, so the sample is of the 2D feature. If εd2 ≪ 1, the coupling between the layers is strong, therefore the sample is of the 3D feature. When εd2 is of the order of 1, the sample shows a layered behavior.
In zero magnetic field, if we take the 2D limit (εd2 ≫ 1) and the 3D limit (εd2 ≪ 1) of Eq. (26) respectively, we can obtain the same expression with the Aslamazov–Larkin theory of the excess conductivity, σs2D ∝ ε−1 and σs3D ∝ ε−1/2.[24,32] The experimental study of dimensionality of conductivity in cuprate superconductor can be found in Refs. [33]–[35].
Figures (a) and (b) show the contour lines of εd2 as a function of temperatures and magnetic fields for BaFe2−xNixAs2 and LaOFeAs. For LaOFeAs, εd2 ∼ 1, it indicates that this material has a layered feature. For BaFe2−xNixAs2, εd2 ≪ 1, we contribute this material to the 3D bulk category.
4. ConclusionWe have investigated the superconducting fluctuation electrical conductivity using SCGT of the Lawrence–Doniach Ginzburg–Landau (LDGL) model. The theoretical results can fit very well with the experimental data of iron-based superconductors in a wide region around Tc. Therefore, the cuprates and iron pnictides superconductors with a layered structure can be well described by the LDGL model. Furthermore, other transport effects of these HTSCs, such as Hall conductivity, the Seebeck effect, and Nernst signal, can also be calculated within the LDGL model.
There are three regions in the H–T phase diagram for cuprates and iron pnictides, i.e., 3D, 2D, and their crossover regions. Although different iron-based superconductors share similar layered structures, the dimensionality could be different. LaOFeAs has strong thermal fluctuations near Tc, so in the H–T phase diagram of this family, the 2D region and the crossover region are relatively large, which proves that LaOFeAs has strong layered or 2D characteristics. BaFe2−xNixAs2 has less fluctuation near Tc, so the 2D region and the crossover region are quite small. BaFe2−xNixAs2 is therefore more of a 3D feature. Considering the dimensionality of LiFeAs discussed in the introduction, we can expect it to show a 3D feature near Tc. For a more specific discussion about the dimensionality of LiFeAs above
, one should use SCGT.